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ABSTRACT 

The spatial emission from the core of cooling flow clusters of galaxies is inadequately de- 
scribed by a (3— model (Cavaliere & Fusco-Femiano 1976). Spectrally, the central region of 
these clusters are well approximated with a two-temperature model, where the inner temper- 
ature represents the multiphase status of the core and the outer temperature is a measure of 
the ambient gas temperature. Following this observational evidence, I extend the use of the 
f3— model to a two-phase gas emission, where the two components coexist within a boundary 
radius r coo i and the ambient gas alone fills the volume shell at radius above r coo i. This simple 
model still provides an analytic expression for the total surface brightness profile: 
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(Note in the first term the different sign with respect to the standard f3— model). Based upon a 
physically meaningful model for the X-ray emission, this formula can be used (i) to improve 
significantly the modeling of the surface brightness profile of cooling flow clusters of galaxies 
when compared to the standard j3— model results, (ii) to constrain properly the physical char- 
acteristics of the intracluster plasma in the outskirts, like, e.g., the ambient gas temperature. 
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1 INTRODUCTION 



To constrain the physical parameters of extended X-ray sources 
(e.g. groups and clusters of galaxies), the observed surface bright- 
ness can be either geometrically deprojected or, more simply, fitted 
with a model obtained from an assumed distribution of the gas den- 
sity. 

Given the hydrostatic equilibrium within the cluster, the grav- 
itational potential supports both the gas and the galaxies distri- 
bution. If the latter is approximated via the King approximation 
(1962) to the inner portions of an isothermal sphere (Lane-Emden 
equation in Binney & Tremaine 1987), the gas density is then writ- 
ten as: 



Assuming isothermality and a /3-model for the gas density 
(eq. [I]), the surface brightness profile has an analytic solution 
(eq. 3.196.2 in Gradshteyn and Ryzhik 1965): 



S(b) = nSr c A(T gas )B(3 / 3-0.5,0.5) 
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Pgas = Po(l + X 



-3/3/2 



(1) 



where x = r jr c and r c is the core radius of the distribution. 

The surface brightness profile observed at the projected radius 
b, S(b), is the projection on the sky of the plasma emissivity, e(r): 



S(b) = 



edr 2 



b 2 \/r 2 — b 2 
The emissivity is equal to 

e(r) = A(T gas ) n v erg s -1 cm -3 , 

where n p = /9 gas / (2.21 /im p ) is the proton density and the cooling 
function, A(T gas ), depends upon the mechanism of the emission 
(mainly due to bremsstrahlung at T gas > 2.5 keV). 



(2) 



(3) 



where the validity of the beta function B(a, b) puts the strict con- 
straint 3/9 > 0.5 and the cooling function A(T gas ) does not change 
radially. 

The fj— model (Cavaliere & Fusco-Femiano 1976, 1978) pro- 
vides a good representation of the observed surface brightness and 
has the advantage to easily constrain the gas density distribution. 

Elsewhere (Ettori 2000), I consider the effect of the presence 
of a temperature gradient in the estimate of the /3— model parame- 
ters. In this paper, I will focus on the deficiency of the /3— model in 
modeling in a satisfactory way the central emission from cooling 
flow clusters of galaxies. 

The cooling flows (e.g. Sarazin 1988, Fabian 1994) result in an 
enhancement of the gas density in the central region due to the high 
cooling efficiency in the cluster core. Recently, there have been 
attempts to model this excess in emission with a generic double 
(3— model, i.e. the sum of two components of surface brightness 
(Ikebe et al. 1996, Xu et al. 1998, Mohr et al. 1999, Reiprich & 
Bohringer 1999). The correlation between the presence of a cool- 
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ing flow and the necessity for a second /3— model is well indicated 
from this figure: Peres et al. (1998) find that 40 per cent of a flux 
limited sample of 55 clusters of galaxies has a deposition rate of 
more than 100 M© yr _1 ; this is the same percentage of the clus- 
ters in the Mohr et al. sample that are better modeled with a double 
f3— model instead of a single one (18 out of 45 clusters). 

The double j3— model as sum of two S(b) in eq. ^, however, 
is not physically meaningful. In fact, a single isothermal tempera- 
ture is usually assumed for both different density components that, 
therefore, are not in equilibrium. Moreover, data with high spatial 
resolution do not show evidence of a second inner core radius. 

I present in this work a simple geometrical and physical model 
of the emission from cooling flow clusters of galaxies. This model 
relies on recent spectral evidence that the cluster plasma can be de- 
scribed as a gas with two phases, one related to the cooling gas and 
the other to the ambient medium. Assuming that the extended intr- 
acluster gas density, n gas , is well described by a j3— model, I show 
in the following section that an analytic expression for S(b) can 
be obtained to describe the surface brightness from cooling flow 
clusters of galaxies. In Section 3, 1 apply this model to real data of 
clusters with or without cooling flows. This model allows to handle 
the emissivity due to each component. I discuss the physical impli- 
cations of this in Section 4. In Section 5, 1 present some concluding 
remarks. 




2 THE TWO-PHASES EMISSION MODEL 



Recent spectral analyses of cooling flow clusters of galaxies (Allen 
et al. 2000, White 2000) have shown how the spectral capabilities 
of the present instruments are unable to resolve all the fine struc- 
tures of a multiphase gas, allowing just a modeling with a two- 
phase component, one that describes the emission from the central 
cooling gas and the other that takes into account the extended emis- 
sion from the ambient medium. 

These observational results provide us with a simple and nat- 
ural model for the total cluster emission: an inner cold phase con- 
fined within r = r CO oi and overlapping the diffuse, ambient gas 
(see Fig.|l]). 

We assume that the two components coexist within r coo i, 
whereas only the ambient plasma fills the cluster volume shell at 
radius above r coo i. The total cluster emissivity is then e(r,T) = 
fcooi + e am b, where 

£ _ f A(Tcooi) n PiCOO i(r) 2 + A(T amb ) n p , amb (r) 2 , r < r coo i 

|_ A(T amb ) n Pjamb (r) 2 , r > r cool ' ' 

from the definition in eq. ||. 

This simple model provides an analytic expression for the sur- 
face brightness profile defined in eq. ^: 
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where the integration limits in S coo i(b) contains the boundary of 
the inner region at r — r coo i. 

Now, I integrate the emissivity along the line of sight. Samb (6) 
is still eq. ^. To integrate S coo i(b), one needs the assumption that 
the only scale parameter of the gas density is the dimension of 
the cooling region, r coo i. Considering that we are in the regime 



Figure 1. Two-phases emission model. r coo i is the boundary of the inner 
region and the only scale parameter used in the fit. 



(r/r coo i) < 1, I can move the sign '-' from the exponent to the 
radix and derive a ft— model in the from of 



Wp.cool ^O.cool 



1 



^cool 



l-5/3 coo i 



(7) 



The behavior of this profile ensures that the gas density within the 
cooling region has no other parameter scale than the dimension of 
the region itself and e coo i goes to zero when r — > r coo i. 

Then, I can integrate analytically S coo \(b) (eq. 3.196.3 in 
Gradshteyn and Ryzhik 1965): 



S COO l(b) — -Scoo^O) 



(— y 

V r coo i / 



0.5+3/5 cool 



Scoot (0) = no,cool'"coolA(Tcool)B(0.5, 3/?cool + 1), 



(8) 



In the equations above, T CO oi and T am b represent the two gas 
temperatures corresponding to the cooling region and to the ambi- 
ent of the cluster, respectively. 



3 COMPARISON WITH THE DATA 

I have applied this model to observations of clusters of galax- 
ies that can map the emission in regions well beyond the clus- 
ter core to disentangle the effect of the two components. More- 
over, I have considered clusters with evidence of a large cooling 
flow and the Coma cluster (ROR: rp800005nOO, exposure time: 
20.0 ksec, 2 = 0.0232), a well-known example of a no-cooling 
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flow cluster. In particular, I have analyzed, as described in Ettori 
and Fabian (1999), the ROSAT PSPC images of A1795 (ROR: 
rp800105n00, exposure time: 33.3 ksec, z = 0.062) and A2199 
(ROR: rp800644n00, exposure time: 33.9 ksec, z = 0.030), that 
have a deposition rate larger than 100 M® yr _1 (Allen et al. 2000) 
and present in the literature detailed spectral analyses of the ASCA 
dataset to be used as reference. 

In Table 1, I quote the results obtained by fitting the az- 
imuthally averaged profiles of the surface brightness between 
and Rout, where the brightness value is larger than 3 times the 
uncertainty in that radial bin. I perform the following fits (see 
Fig. ^): (i) single /3— model, (ii) the double (3— model presented 
here with 6 parameters [S coo i(0), r coo i, /3 coo i, S am b(0), r c , /3 amb ], 
(iii) the double f3— model with 5 parameters, i.e. fixing r c = r coo i. 

The decrease of the \ 2 is significant where a cooling flow is 
present. Where this decrease is not meaningful (or not present), 
like in Coma cluster, the single j3— model still represents a good 
description of the data. 

Where a two-phase model provides a significantly better fit, 
I find that the F-test shows no statistical improvement with a 6 
parameters fit (cf. Table 1). Therefore, I use in the following con- 
siderations the best fit results obtained with a 5 parameters fit, i.e. 
r c = r coo \. This is not in contradiction with the present observa- 
tional results. Allen (2000) quotes the cooling radii obtained from 
deprojection analysis of 30 cooling flow clusters images. Ettori & 
Fabian (1999) estimate the core radii for 23 of these clusters using 
a single j3— model over the radial range [0.1, 1] rsoo, the radius at 
which the mean cluster density is 500 the background value. The 
distribution of the ratio, r c /r coo \, has a median value of 1.33, an 
average of 1.61 and a standard deviation of 1.19, and can be con- 
sidered consistent with ~ 1. For the clusters in exam here, I mea- 
sure a r c /r coo i ratio of 1.32(±0.24, 90 per cent confidence level) 
and 0.84 (±0.15) for A1795 and A2199, respectively. I remind, 
however, that I am using a different definition of r CO oi than the one 
adopted in the standard spatial analysis: in the latter, r coo i is the 
radius where the cluster cooling time first exceeds the Hubble time, 
whereas in this work r CO oi defines the boundary of the central cool 
phase of the gas. 

It is worth to note that another version of a '5 parameters' fit, 
in which /3 coo i is fixed equal to /3 am b and r c , r coo \ are left free to 
vary, provides a significantly worse x 2 man me one obtained by 
using the '5 parameters' fit adopted here. 

Finally, the \ 2 obtained with the models above does not vary 
significantly from the % 2 measured after the fit with other mod- 
els which are not strictly based upon a physical framework, like 
the sum of two standard j3— models with 6 free parameters (e.g. 
Reiprich & Bohringer 1999) or, fixed the slope to a common value, 
with 5 free parameters (Mohr et al. 1999). 



n p (Q) 



4 DISCUSSION 

The use of this physically meaningful model allows us to directly 
handle each of the two gas distributions, one that describes the pro- 
file of the gas related to the cooling flow and the other that model 
the ambient gas. As shown above, there is no statistical justifica- 
tion for using the 6 parameter fit. Therefore, I consider hereafter 
the case r c — r coo \. 

We use in the following discussion the definition of the central 
gas density for a cluster at redshift z, that is 



4tt(1 + z) 4 S(0) 



1.21r c A(T)5 



(9) 



where %>(0) is in cm -3 , r c is the core or cooling radius in cm, 
S(0) is the central surface brightness in cts s" 1 ster -1 , 5 is the 
proper beta function and the cooling function A(T) is in unit of cts 
s _1 cm 5 . 

Using the condition that the two phases have to be in pres- 
sure equilibrium, I can now put constraints on their temperatures. 
To do this, I consider the mean properties of each phase to handle 
integrated values instead of differential ones, because of the simple 
assumption that each phase is represented by a single temperature 
that does not depend on the radius. Therefore, each phase density, 
n(< r), is the integral of the radial density, n(r), over the volume 
occupied from that phase (between and r coo i for the inner phase; 
between and r — Xr coo \ for the outer component) divided by the 
integrated volume. Then, 



T coo \ — Xai 



^amb 
T^cool 



= Iamb g 



l/(2-c0 r2/(2-a) 



= T amb /(10) 



where I have made use of the relation in eq. ^ A(T)riQ ~ 
S(0)/B(a,b), I have assumed A(T) ~ T a (a w 0.5 for only 
bremsstrahlung emission observed by broad-band instruments), 
and I have defined 

A(T am b) ^ amb 

(0) _ 5 am b(0) Bcooi 



9 = 



A(T C ooi)n coo i(0) 2 S'cooi(O) B amb 
?"t am b(< r) n coo i(0) 
n am b(0) n coo i(< r CO oi) 

1 j*(l + x 2 )- 1 ' 5,3 ^x 2 dx 

x3 J^(i-x 2 y- 5 ^x 2 dx 

Ix[l-5]/X 3 



0.5B(1.5, 1.5/3 C ooi + 1) 

/ = ff V(2"«) J 2/(2-«) j (U) 

with /^(l - x 2 ) 1 - 5/3 -°° l x 2 dx = 5(1.5, 1.5/3 C ooi + l)/2 and 

Ix[a] = f X (l + x 2 )- a ^x 2 dx. 

However, one generally measures a single emission-weighted 
temperature, T. Given the considerations above, I can now disen- 
tangle the two components (if any) using eq. |nj in the following 



relation: 



T 



= T; 



1+/ B(1.5, 3&00I + I) (2gl x [3])- 1 



-5(1.5,3/30001 + 1) (2gIx[S])- 



^amb F. 



(12) 



where I still use the relation A(T)rio ~ S(0)/B(a,b), calculate 
^(l-x 2 ) 3l3 -°° l x 2 dx = 5(1.5, 3/3 coo i + l)/2 and adopt the sym- 
bols f,g,Ix defined in eq. [ill. 

In the equations above, both the function / and F have to 
be smaller than 1 by definition. Their behaviour, however, depends 
strongly upon X, the radius in unit of r coo i up to where the outer 
phase extends and can be represented with a single temperature. 
Figure |3| shows how the function / and F depend upon X: f dimin- 
ishes significantly due to the presence of X 3 in I (eq. pi] ), whereas 
F converges quite rapidly (at X J> 4), providing a robust estimate 
on the T/T am b ratio. Therefore, even if we are not able to constrain 
the ratio between the temperature of the two phases due to the un- 
certainty of the extension of the outer component, we can assess the 
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Table 1. Best-fit results using the models discussed in the text. The core (cooling) radii are in h 50 Mpc. The values for the F-test represent the level of 
significance. The symbols 20% and 20a indicate the two 0— model with 5 and 6 parameters, respectively. 



cluster Rout 
Mpc/' 



10 







2 0-5 params 
)"cooi,X 2 (d-o.f.) 

/^cool i ^amb 



2/3-6 params 
f C ool,^c,X 2 (d-o.f.) 

/^cooli ^amb 



10 



F-test 
■ 2/3 5 2/3 5 



20e 



A1795 1.49/15.2 



A2199 1.50/30.2 



0.10 ±0.01, 40.9(28) 

0.631 ±0.002 
0.08 ±0.01, 23.1(58) 
0.586 ±0.001 
Coma 1.49/38.2 0.47 ± 0.01, 8.2(74) 
0.801 ± 0.008 



0.26 ±0.01, 8.6(26) 
1.801 ±0.126, 0.761 ±0.009 

0.13 ±0.01, 4.7(56) 
1.635 ± 0.058, 0.644 ± 0.004 

0.51 ±0.01, 6.5(72) 
9.071 ± 1.532, 0.835 ± 0.009 



0.64 ± 0.15, 0.28 ± 0.01, 8.2(25) 

4.459 ± 1.012, 0.780 ±0.013 > 0.99 0.45 
0.74 ± 0.06, 0.13 ± 0.01, 3.6(55) 

67.2 ±9.9, 0.650 ±0.005 > 0.99 0.85 
33.67 ± 3.05, 0.60 ± 0.03, 8.6(71) 

572.4 ± 109.1, 1.376 ±0.116 0.83 0.12 




Figure 2. The cluster surface brightness profile is here fitted with a single 0— model (dashed line) and a 2 0— model (solid line; the thickest indicates the 
model with 5 parameters). The dotted lines show the two components of the 2 0— model with 5 parameters. The panel below shows the residuals in unit of a. 



ambient temperature, T amD , in a cooling flow cluster with an emis- 
sion weighted temperature, T, just using the azimuthally averaged 
surface brightness profile. 

I show in Table 2 the constraints on F = T/T am b obtained 
from the spatial fit of the surface brightness profiles of the clusters 
in exam and compare these values to the two-temperature spectral 
results in Allen et al. (2000), Markevitch et al. (1998, 1999) and 
White (2000). The agreement is remarkably good with the results 
of Markevitch and collaborators and White, which assume a two- 
phases gas for their spectral model in a way similar to the one I 
have adopted for the physical framework described above. (Note 
that Markevitch et al. measure the isothermal ambient temperature 
excising the cooling region, whereas White adds a cooling flow 
component in the spectral fit). On the other hand, the disagreement 
with the results of Allen and collaborators can be explained with the 
more complex model that they adopt, where an absorption intrinsic 
to the cluster is combined with the cooling flow only. 

Several aspects of the cluster physical characteristics are af- 
fected from the inclusion of a cooling flow in the modeling of 
the surface brightness with a /3— model. With respect to the single 
/3— model, one expects (i) excess in the gas density in the cooling 



region, (ii) change in the /3 amD value, (iii) variations in the gas am- 
bient temperature. 

I present in Table 3 some of the more interesting physical 
quantities that can be evaluated with the equations above and given 
an emission-weighted temperature, T (from Allen et al. 2000: 
kT = 5.40 ± 0.05 and 4.16 ± 0.03 keV for A1795 and A2199, 
respectively). 

For example, if one identifies the inner component with the 
cooling flow, a proper description of its gas distribution is now 
available through eq. ^. The luminosity of the intracluster gas can 
then be estimated without the contribution of the emission from the 
cool phase: 

Lamb = J 1.21np(r) A(T amb ) 4nr 2 dr (13) 

where the integral is computed upon the cluster volume and A(T) 
is here in erg s -1 cm 3 . Using only the cluster surface brightness 
profile and a broad-band emission-weighted gas temperature, and 
applying eq. |l|and|l| I will investigate in a forthcoming paper the 
effects on the clusters luminosity-temperature relation of the pres- 



© 0000 RAS, MNRAS 000, 000-000 



5 



Table 2. Comparison between the ratio of the gas temperature obtained spectrally from (i) a single temperature model and (ii) a two-temperature (or isothermal 
+ cooling flow components) model and the F function in eq. [l^. The errors are at the 68.3 per cent confidence level. In the observational data, these come from 
the propagation of the error on the temperature estimates. The error on the F function is computed after 500 Monte Carlo fit of the original surface brightness 
profile. References for the observational data: [1] Allen et al. (2000), [2] Markevitch et al. (1998, 1999), [3] White (2000). 



cluster [1], [2], [3] F: best-fit results 



A1795 0.870 ± 0.023, 0.769 ± 0.065, 0.799 ± 0.057 0.726 ±0.010 
A2199 0.837 ±0.037, 0.917 ±0.035, 0.909 ±0.030 0.905 ± 0.006 



Table 3. Ratios between the quantities discussed in the text and estimated by using (i) a fit with a 1 — fi and T and (ii) a 2 — /3s model and T amb , respectively. 
The cooling function is estimated using a MEKAL model (Kaastra 1992, Liedhal et al. 1995) in XSPEC (version 10, Amaud 1996). The typical relative errors 
on the ratio of luminosities, gas and total masses are of about 1 per cent and come from 500 Monte Carlo replications of the observed surface brightness profile 
and emission-weighted temperature. 



cluster A(T)/A(T amb ) L/L amb M g as/M gas Mtot/AAot 

0.5, 1.0, 1.5 Mpc 0.5, 1.0, 1.5 Mpc 0.5, 1.0, 1.5 Mpc 



A1795 0.87 ±0.01 1.34,1.18,1.16 0.89,0.92,0.97 0.73,0.64,0.62 
A2199 0.96 ±0.01 1.06,1.04,1.05 0.94,0.98,1.02 0.86,0.84,0.83 



1.4 
1.2 
1.0 
0.8 

<4-T 

0.6 

0.4 

0.2 
0.0 







!C3. 

k 



□ 



□ 



□ 
-□ 



10 



12 



X (r/r oool ) 



Figure 3. Plot of the functions / (diamonds) and F (squares) vs. X. The 
dashed lines represent A1795, the dotted ones A2199. The error-bars are 
overplotted to each point. 



where A/ gas , coo i = 2ty p , coo i rj? ool 73(1.5, 1.5/3 coo i + 1) w 
27 (3) x 10 n M© for A1795 (A2199) for an assumed T coo i = 
(f/F) T = 0.2T; the total gravitating mass is given by the appli- 
cation of the hydrostatic equilibrium, 



M t , 



T^amb Pamb ) 



fj,m p Gp g!lB 



dr 



(15) 



and is proportional to r coo i/3 amb Tamba; 3 (l 



at r > r coo i. 



5 CONCLUSIONS 

I have presented a new analytic formula to model the total surface 
brightness profile of clusters of galaxies where a two-phases intr- 
acluster gas can be assumed. This scenario is consistent with the 
present results of spectral analyses of the central regions of clusters 
that harbour a cooling flow. 

The use of this formula allows to properly disentangle the con- 
tribution of the cooling flow to the cluster emissivity using only 
the spatial distribution of the X-ray photons. After removing the 
contamination from the cooling flow, I show how some relevant 
physical parameters are affected, like, for example, the ambient gas 
temperature (see Table 2). In a forthcoming paper, I will investigate 
the systematic changes in the temperature, luminosity and mass (cf. 
Table 3) of a sample of clusters of galaxies and how these variations 
affect the relations among these quantities. 



ence of significant cooling flows (see, e.g., the results from spectral 
analyses in Allen & Fabian 1998, Markevitch 1998). 

Appreciable corrections can also affect M gas , M to t and the 
terms of the so-called f3— problem (Mushotzky 1984, Edge & Stew- 
art 1991) due to the variation of the (5 value (for A1795 and A2199, 
/3amb increases by 20 and 10 per cent when compared to the 1 — /3 
model fit results, respectively). In the two-phases model described 
here, it is simple to calculate the gas and the total mass: M gas is the 
integral of p gas = p coo i + p am b = 2.21^m p (n p , coo i + n Piam b) 
upon the cluster volume and, in particular, 



M gas (r > r coo i) = M gas ,cooi + / Pamb 47rr 2 dr 



(14) 
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